Cosmology with Hu-Sawicki gravity in Palatini Formalism 
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Cosmological models based on /(i?)-gravity may exhibit a natural acceleration mechanism without 
introducing a dark energy component. In this paper, we investigate cosmological consequences of 
the so-called Hu-Sawicki f{R) gravity in the Palatini formalism. We derive theoretical constraints 
on the model parameters and perform a statistical analysis to determine the parametric space 
allowed by current observational data. We find that this class of models is indistinguishable from 
the standard ACDM model at the background level. Differently, from previous results in the metric 
approach, we show that these scenarios are able to produce the sequence of radiation-dominated, 
matter-dominated, and accelerating periods without need of dark energy. 
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I. INTRODUCTION 

General Relativity (GR) is a very well- tested and es- 
tablished theory of gravity. However, all the successful 
tests performed so far do not account for the ultra-large 
scales corresponding to the low curvature characteris- 
tics of the Hubble radius today. Thus, in principle it 
is conceivable that the recent discovery of the late-time 
cosmic acceleration [l| could be explained by modifying 
Einstein's GR in the far infrared regime. Although the 
diversity of approaches is the hallmark in this field (see, 
e.g., Refs. [2]), the simplest possible theory results by 
adding terms proportional to powers of the Ricci scalar 
R to the Einstein-Hilbert Lagrangian, the so-called f{R) 
gravity (see Refs. Q for recent reviews). 

Differently from the standard general relativistic sce- 
nario, f{R) cosmology can naturally drive an accelerat- 
ing cosmic expansion without introducing dark energy. 
However, the freedom in the choice of different func- 
tional forms of f{R) gives rise to the problem of how 
to constrain the many possible f{R) gravity theories. 
In this regard, much efforts have been developed so far, 
mainly from the theoretical viewpoint [4]. General prin- 
ciples such as the so-called energy conditions [5], non- 
local causal structure [6], and issues such as loop quan- 
tum cosmology |7|, have also been taken into account 
in order to clarify its subtleties. More recently, obser- 
vational constraints from several cosmological data sets 
have also been explored for testing the viability of these 
theories 

An important aspect that is worth emphasizing con- 
cerns the two different variational approaches that may 
be followed when one works with f{R) gravity, namely, 
the metric and the Palatini formalisms. In the metric 
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formalism the connections are defined a priori as the 
Christoffel symbols of the metric and the variation of 
the action is taken with respect to the metric, whereas 
in the Palatini variational approach the metric and the 
affine connections are treated as independent fields and 
the variation is taken with respect to both (for a review 
on f{R) theories in the Palatini approach see [9]). The 
result is that in the Palatini approach the connections de- 
pend on the metric and also on the particular f{R), while 
in the metric formalism the connections depend only on 
the metric. We have then that the same f{R) leads to 
different spacetime structures. 

These differences also extend to the observational as- 
pects. For instance, we note that some cosmological mod- 
els based on a power-law functional form in the met- 
ric formulation fail in reproducing the standard matter- 
dominated era followed by an acceleration phase 
(see, however, [TH), whereas in the Palatini approach, 
some analysis have shown that such theories admit the 
three post inflationary phases of the standard cosmolog- 
ical model [12]. Nevertheless, we do not have yet a clear 
comprehension of the properties of the Palatini formula- 
tion of f{R) gravity in other scenarios, and issues such 
as the Newtonian limit [13] and the Cauchy problem 
are still contentious (see, however, Refs. [15|). Another 
issue with the Palatini f{R) formulation concerns surface 
singularities of static spherically symmetric objects with 
polytropic equation of state [9|, ll6|. This problem has 
been reexamined in Refs. [l^, where the authors showed 
that the singularities may have more to do with pecu- 
liarities of the polytropic equation of state used (e.g. its 
natural regime of validity) and the f{R) model chosen. 

Although being mathematically simpler than the met- 
ric formulation, the Palatini approach has received little 
attention from the point of view of cosmological tests. 
Following early works in this direction [l2|, [ill, [l9|, in 
this paper we explore the cosmological consequences of 
a class of modified f{R) gravity, recently proposed by 
Hu & Sawicki [20|, in the Palatini formulation. Among a 
number of f{R) models discussed in the literature, this 
is designed to posses a chameleon mechanism that allows 
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to evade solar system constraints and the cosmological 
scenario that arises from the metric formahsm of this 
model has been shown to satisfy the conditions needed 
to produce a cosmologically viable expansion history. In 
order to test the observational viability of this class of 
models in the Palatini approach, we use different kinds 
of observational data, namely: type la supernovae (SNe 
la) observations from Union2.1 sample [21|, estimates of 
the expansion rate at z ^ 0, as discussed in Ref. [22], 
current measurements of the product of the CMB acous- 
tic scale and the comoving sound horizon at photon 
decoupling (CMB/BAO ratio) [23] and measurements of 
the gas mass fraction in galaxy clusters [24]. We show 
that for a subsample of model parameters, the so-called 
Hu-Sawicki scenarios in the Palatini approach are indis- 
tinguishable from the standard ACDM model, being able 
to produce the sequence of radiation-dominated, matter- 
dominated, and accelerating periods without need of dark 
energy. 



pressure p whose energy-momentum tensor is given by 
Tt^ = diag{-p,p,p,p). 

The generalized Friedmann equation, obtained from 
(|2]) and (|3]), can be written in terms of the redshift z as 



^0 
where 



2 fn RJrr — fn 



, (4) 



(5) 



and ^rno and QrO stand for the present-day values of the 
matter and radiation density parameters. The trace of 
([2al) results in 



RfR - 2/ = -3H^nmo{l + zf , 



(6) 



which we will consider in order to restrict the parameters 
of the given f{R) theory. 



II. THE PALATINI APPROACH 



The modified Einstein-Hilbert action that defines an 
f{R) gravity is given by 



S- 



-gf{R) + 



(1) 



where hc'^ = SttG, g is the determinant of the metric ten- 
sor and Sm is the standard action for the matter fields. 

In the Palatini variational approach the fundamental 
idea is to regard the torsion- free connection Fj^^, as well 
as the metric g^i, entering the action above, as two inde- 
pendent fields. The field equations in this approach are 



and 



V/9 (v^/i?^^^) = 0, 



(2a) 



(2b) 



where fn = df/dR and V represents covariant deriva- 
tives. From the above equations, we obtain the connec- 
tions 



r)fR, (3) 



where {^j^} are the Christoffel symbols of the metric. As 
mentioned earlier, the connections, and so the gravita- 
tional fields, are described not only by the metric ^^j, 
but also by the proposed f{R) theory. As a cosmologi- 
cal model we consider a homogeneous isotropic Universe 
described by the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) flat geometry ^^j, = diag{—l^o?^o?^o?)^ where 
a{t) is the cosmological scale factor, and as the source 
of curvature a perfect-fluid with energy density p and 



A. The Hu-Sawicki gravity model 

In this paper we are interested in the f{K) modified 
gravity model proposed by Hu and Sawicki [20| 



f{R) = R-m^ 



^1 {^r 



C2{ 



R 



(7) 



For the case n = 1, the above equation can be rewritten 

as m 



f{R) = R 



2A 



where A = m?ci/2c2 and pi? = m? /c2 are dimensionless 
parameters and m? = H^VtrnO' Note that, in the regime 
R » /i^, the model above is practically indistinguish- 
able from the ACDM scenario^. We have verified that the 
parameter n is completely unconstrained by current data. 
Thus, without loss of generality and following Ref. [26|, 
we focus our analyses on the case n = 1. 

The model above is intended to explain the current 
acceleration of the universe without considering a cos- 
mological constant term. Cosmological and astrophysi- 
cal constraints on the Hu-Sawicki gravity scenario in the 
metric formalism have been examined in a number of pa- 
pers (see, e.g., Refs. \27]). However, there is a complete 
lack of studies on the mechanism of this proposal in the 
Palatini approach. Note that, irrespective of the formal- 
ism adopted, the model ([7]) must satisfy certain viability 



^ As shown in Ref. [26], this is equivalent to have |/ho| ~ 1? where 
the subscript denotes present-day quantities. See also Table I 
for some estimates of this quantity. 
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FIG. 1: (left) The predicted Hubble evolution H{z) as a function of the redshift for the Hu-Sawicki model in the Palatini 
formalism (Eqs. [4][7|. The curves correspond to the best-fit values of Qmo and C2 discussed in the text. For the sake of 
comparison, the standard ACDM model prediction is also shown. The data points are the measurements of the H{z) given in 
Ref. [23] . (right) Hubble diagram for 580 SNe la from the Union2.1 sample 21]. The curves correspond to the best-fit values 
of Qmo and C2 displayed in Table [H 



conditions, as for example, the positivity of the effec- 
tive gravitational coupling which requires / fn > (to 
avoid anti-gravity). For the current epoch, this condition 
implies 



3 + 2(^)2c2 



(8) 



which corresponds to the following bounds for the pa- 
rameter 02'. 



III. OBSERVATIONAL CONSTRAINTS 

In order to perform the observational tests discussed 
below, we solve Eq. (|4]) by taking the following steps: i) 
we first set z = and compute ci (in terms of C2 and Rq) 
from Eq. ([6]); ii) we combine the latter result with Eq. (j4]) 
in order to obtain a solution for Rq] iii) the value of Rq is 
then used as initial condition to solve Eq. ([6j) by using a 
fourth-order Runge-Kutta method and, finally, to obtain 
the function H{z)^ as given by (j4]). 

Fig. [1] shows the evolution of the Hubble param- 
eter and the predicted distance modulus jj.{z) = 
51og[(iL(^)/Mpc]+25, where dL = (1 + ^) Jq ^{z') stands 
for the luminosity distance, as a function of redshift for 
some best-fit values for VtrnQ and C2 obtained in our anal- 
yses (see Sec. IIIC). For the sake of comparison, the stan- 
dard ACDM prediction with VtrnQ = 0-26 is also shown 
(dashed line). 



A. Expansion rate 

In order to test the observational viability of the f{R) 
scenario described by Eq. ([7j) , we firstly use current mea- 
surements of the expansion rate at z 7^ 0. Currently, 
this test is based on the fact that Luminous Red Galax- 
ies (LRGs) can provide us with direct measurements of 
H{z) [28] (see also [29| for a recent review on H{z) mea- 
surements from different techniques). This can be done 
by calculating the derivative of cosmic time with respect 
to redshift, i.e., H{z) 0^ ~(T+i)'^- This method was 
first presented in [28| and consists in measuring the age 
difference between two red galaxies at different redshifts 
in order to obtain the rate Az/At. By using a recent 
released sample from Gemini Deep Survey (GDDS) (3Q| 
and archival data [31], Ref. [22] has calculated 11 H{z) 
data points over the redshift range 0.1 < z < 1.75. 

We estimate the free parameters P by using a statis- 
tics, i.e., 



Xh 



Nh 
i=l 



'Hk{zi\p)-Hi,^{z,)y 



(10) 



where i^*^(z|P) is the theoretical Hubble parameter at 
redshift Zi and ai is the uncertainty for each of the 
Nh = 11 determinations of H{z) given in [22]. In our 
analyses, the Hubble constant Hq is considered as a nui- 
sance parameter and we marginalize over it. 



B. SNe la 

The predicted distance modulus for a supernova at 
redshift z, given a set of parameters P, is /i^(2:^|P) = 
5 log 1^1,(^2, P)+/^o 5 where Dl = HqcIl is the Hubble- free 
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FIG. 2: (left) Contours of Ax^ — 2.3 and Ax^ = 6.17 in the Qmo x C2 plane arising from SNe la union2.1 sample [2H 
(Orange) and current H{z) data [l^] (Green). Black contours stand for the joint analysis involving these two sets of data. 
{right) Contours of Ax^ = 6.17 in the Q^mo x C2 plane when 7 measurements of the CMB/BAO ratio and 57 measurements of 
the gas mass fraction of galaxy clusters are added. We have also indicated the region Jrq < corresponding to the bounds JS} 
for the best-fit values shown in Table I. 



TABLE I: Best-fit values for C2, ci and /rq. The error bars for Qmo correspond to 68.3% and 95.4% confidence intervals. 



Test 


QrnO 






C2 


Cl 


Iro 


\min 


Xr 


H{z) 


0.22 


+0.11 

-0.08 


+o.iy 

-0.13 


0.08 


2.23 


0.88 


7.72 


0.86 


SNe la 


0.29 


+0.19 
-0.17 


+0.34 
-0.21 


-0.40 


-5.41 


1.03 


553.81 


0.96 


CMB/BAO 


0.29 


+0.03 
-0.01 


+0.07 
-0.03 


-0.18 


-2.17 


1.08 


1.24 


0.25 


fgas 


0.28 


+0.01 
-0.01 


+0.04 
-0.02 


-0.03 


-0.04 


1.61 


95.10 


1.73 


H{z) + SNe la 


0.28 


+ 0.06 
-0.07 


+ 0.12 
-0.15 


-0.85 


-12.64 


1.02 


561.63 


0.95 


CMB/BAO + fgas + H{z) + SNe la 


0.30 


+ 0.01 
-0.01 


+ 0.02 
-0.02 


-0.25 


-3.04 


1.06 


659.87 


1.01 



luminosity distance, fio = 42.38 — Slog^^g ^ ^ stands 
for the Hubble constant Hq in units of 100 km/s/Mpc. In 
our analysis, we use the latest Union2.1 sample which in- 
cludes 580 SNe la over the redshift range 0.1 < z < 1.75. 
In order to include the effect of systematic errors into 
our analyses, we follow [21] and estimated the best fit 
and error bars to the set of parameters P by using a 
statistics, with 

580 

Xlive = E(/^°'' - MD(C-^),,(Mf - , (11) 

where C is a 580 x 580 covariance matrix [21]. As for the 
analysis involving H{z) data, we also marginalize over 
Ho. 

C. Results 

In Fig. [2] (left) we show the first results of our sta- 
tistical analyses. Contours corresponding to Ax^ =2.3 
and Ax^ = 6.17 in the flrno x C2 plane are shown for the 
given above. From this analysis, we clearly see that 
neither the current H{z) and SNe la measurements alone 



nor the combination of them (x^ = XsNe~^XH) can place 
tight constraints on the values of the f{R) parameter 
C2 . We also observed that other cosmological observables 
seem to be unable to provide orthogonal contours to those 
shown in Fig. [2] (left), which in principle would lead to 
tighter bounds on C2 from a joint analysis. Fig. [2] (right) 
shows the results obtained when two other sets of cosmo- 
logical observations are included, namely: 7 estimates of 
the ratio of the CMB acoustic scale Ia and the baryonic 
acoustic oscillation (BAO) peak, as given in Ref. p3| , 
and the 57 measurements of the gas mass fraction in X- 
ray luminous galaxy clusters discussed in Ref. Q (we 
refer the reader to Refs. |25| for more details on these 
cosmologial tests). For this joint analysis, the best-fit 
values are ftrno =0.30 and C2 = —0.25, with the reduced 
Xr = Xmin/^ — ^'^^ defined as degrees of freedom), 
which is very close to the value found for the standard 
ACDM model, Xr — 1-001. This result clearly shows that 
these scenarios are indistinguishable from each other at 
the background level, which coincides with the analy- 
sis performed in the metric formalism and discussed in 
Ref. [34] (although the expansion histories for both met- 
ric and Palatini formalisms are completely different). It is 
worth mentioning that in our analysis we consider both 
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FIG. 3: Deceleration parameter (left) and effective equation of state (right) as a function of z for the best-fit values of Qmo 
and C2 presented in Table [H The ACDM model is also shown for the sake of comparison. 



positive and negative values of the f{R) parameter C2, 
in agreement with the constraint ([9]). Note also that for 
C2 = 0, the cosmological scenario discussed above reduces 
to the Einstein-de Sitter model, being ruled out by cur- 
rent data. We use the joint best-fit value discussed above, 
as well as the others shown in Table HI to discuss some 
cosmological consequences of this class of f{R) cosmolo- 
gies in the next section. 

IV. COSMOLOGICAL CONSEQUENCES 

Two important quantities directly related to the ex- 
pansion rate of the Universe and its derivatives are the 
deceleration parameter 

q{z) = ^^^H'{z) - 1 , (12a) 

and the effective equation of state (EoS) 

«'e// = -l + ^^^if'(^), (12b) 

where a prime denotes differentiation with respect to z 
and H{z) is given by Eq. (j4]). For all sets of best-fit val- 
ues displayed in Table HI we show q{z) and Weff{z) as a 
function of the redshift in Fig. [3l Similarly to the analy- 
ses presented in Sec. Ill, we have included a component 
of radiation {ftro = 5 x 10~^) to plot these curves. As 
can be seen from Fig. [3] (left), for some combinations of 
parameters the cosmic evolution is well-behaved, with a 
past cosmic deceleration and a current accelerating phase 
starting at z :^ 1, which seems to be in good agreement 
with current kinematic analyses (see, e.g., [33]). 

For all best-fit combinations shown in Table HI 
Fig. [3] (right) shows that the Universe went through 
the last three phases of cosmological evolution, i.e., 
a radiation-dominated {weff = 1/3), a matter- 
dominated {weff = 0) and a late-time accelerating 



phase {weff < 0), which is similar to what happens in 
the standard ACDM cosmology. From these results, it 
is also clear that the arguments of Ref. [10] about the 
behavior of Wefj in the metric approach (we refer the 
reader to pjj] for a different conclusion) seems not to 
apply to the Palatini formalism, at least for the class 
of models discussed in this paper and the interval of 
parameters ftmo and C2 given by our statistical analysis 
- the same also happens to power-law [l2|, IHI and 
exponential [19] f{R) gravities in the Palatini formalism. 
Note also that, differently from exponential- type f{R) 
models (see, e.g., [HI), we have not found solutions of 
transient cosmic acceleration in which the large-scale 
modification of gravity will drive the Universe to a new 
matter-dominated era in the future. 



V. CONCLUSIONS 

/(i?)-gravity provides an alternative way to explain 
the current cosmic acceleration with no need of invok- 
ing either the existence of an extra spatial dimension or 
an exotic component of dark energy. Among a number 
of f{R) models discussed in the literature, the so-called 
Hu-Sawicki scenarios are designed to posses a chameleon 
mechanism that allows to evade solar system constraints. 
Although the cosmological scenario that arises from the 
metric formalism of this model has been shown to satisfy 
the conditions needed to produce a cosmologically viable 
expansion history, its Palatini version had not yet been 
investigated. 

In this paper, we have discussed several cosmologi- 
cal consequences of this class of models in the Palatini 
formalism. We have performed consistency checks and 
tested the observational viability of these scenarios by us- 
ing one of the latest SNe la data, the so-called Union2.1 
sample, measurements of the expansion rate H{z) at in- 
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termediary and high-z, estimates of the CMB/BAO ratio 
and current observations of the gas mass fraction in X-ray 
luminous galaxy clusters. We have found a good agree- 
ment between these observations and the theoretical pre- 
dictions of the model, with the reduced Xmin/^ — ^ 
the analyses performed (see Table I). In what concerns 
the past cosmic evolution predicted by this class of mod- 
els, we have shown that for a subsample of model param- 
eters, the so-called Hu-Sawicki scenarios in the Palatini 



approach are indistinguishable from the standard ACDM 
model, being able to produce the sequence of radiation- 
dominated, matter-dominated, and accelerating periods 
without need of dark energy. 
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